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Photogalvanic Effects in Heteropolar Nanotubes 
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We show that an electrical shift current is generated when electrons are photoexcited from the 
valence to conduction bands on a BN nanotube. This photocurrent follows the light pulse envelope 
and its symmetry is controlled by the atomic structure of the nanotube. We find that the shift current 
has an intrinsic quantum mechanical signature in which the chiral index of the tube determines the 
direction of the current along the tube axis. We identify the discrete lattice effects in the tangent 
plane of the tube that lead to an azimuthal component of the shift current. The nanotube shift 
current can lead to ultrafast opto-electronic and opto-mechanical applications. 
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Recent progress in synthesis of nanometer scale mate- 
rials has led to the discovery of BxCyN^ nanotubes. They 
can be formed in nearly homogeneous ordered structures 
[Q,^ or in multiwall hybrid structures with alternating 
BN and C compositions |^ . These materials differ funda- 
mentally from the structurally similar carbon nanotubes 
by being non- centra symmetric (NCS) and polar. This 
opens the possibility for accessing a new class of photo- 
voltaic effects at the molecular scale. 

Photovoltaic effects in NCS materials are often based 
on an asymmetric generation of hot carriers at mo- 
menta ±k leading to a so called "ballistic photocurrent" 
[Q-^. In polar NCS materials, photoexcitation across the 
bandgap with polarized light also produces the so called 
electrical "shift current" . Microscopically this orig- 
inates from a net displacement of charge in the unit cell 
due to light induced interband transitions (Je), intraband 
relaxation (Jg) and (radiative) transitions to the original 
bands (Jr)- In bulk materials, the excitation component 
Je usually prevails. 

We show that unpolarized light can induce a shift cur- 
rent Je in polar NCS nanotubes, with a direction along 
the tube axis which is determined by the chiral index 
of the tube. We find that this effect has an essentially 
quantum mechanical origin, where the sign of the current 
along the tube axis is controlled by the phase matching 
of electronic Bloch waves around the tube circumference. 
The discrete lattice structure controls the azimuthal com- 
ponent of this current and thus produces a net helical 
current on the tube. These photo-effects can lead to 
an assortment of new opto-electronic, opto-mechanical 
and magnetic applications, and interesting extensions to 
ring structures ||l^ and heterojunctions ||l^. We should 
stress that the photoeffccts we discuss here are found at 
nonlinear order in the exciting fields, and are therefore 
physically distinct from chiral currents tilted with respect 
to the tube axis, predicted in dc bias driven chiral BN 
nanotubes |^2| or chiral stretched C nanotubes 14|. 



A flat BN sheet has a honeycomb lattice with the B and 
N occupying alternate sublattices as shown in Figure 1. 
The physics of the shift current can be understood quali- 
tatively by considering the response of this 2D network of 
bonds to normally incident polarized light. The valence 
states are polarized towards the N sites and the conduc- 
tion states towards the B sites. For vertically polarized 
(y-polarized) incident light the response of the system is 
dominated by the bonds which produce a net y polarized 
current. For horizontally (x-polarized) incident light, the 
bonds with a nonzero horizontal component dominate, 
but these also produce a net y-polarized current. Note 
that excitation with an unpolarized incident field cannot 
produce a net current on this lattice since it has a three- 
fold symmetry. However, this symmetry is removed when 
the sheet is wrapped into a cylindrical nanotube where 
the physics should be dominated by excitations with the 
field polarized along the tube axis. Depending on the 
wrapped structure of the tube, we anticipate an intrinsic 
photocurrent which can flow along the tube, around the 
tube, or in a chiral pattern on the tube surface. 

These ideas can be quantified by developing a quantum 
mechanical model which generalizes our long wavelength 
theory for a nonpolar carbon nanotube [ p^ to the het- 
eropolar lattice. We work in a basis of Bloch orbitals 

where a = ±1 denote the two sublattices occupying sites 
at Tq in the n-ih unit cell. We study the states near 
the conduction and valence band edges at the K and K' 
points of the Brillouin zone shown in Figure 1. A long 
wavelength Hamiltonian for these states is obtained by an 
expansion in small crystal momenta around these points, 
and yields in our two component basis (fi = 1) 



A Xvpik — iSx) 

XvF{k + iSx) -A 



(1) 



where A = ±1 is an index which labels an expansion 
around the K or K' points. The Hamiltonian in Eq. (|^) 
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is parameterized by three energies: a symmetry breaking 
site diagonal potential A, the kinetic energy vpk for mo- 
tion along the tube, and the kinetic energy vpS due to 
the quantization of the transverse momentum around the 
circumference of the tube. On a cylindrical tube with lat- 
tice constant a and primitive lattice translation vectors 
in its tangent plane Ti — a(l,0) and T2 — a(l/2, \/3/2) 
(Figure |l|) it is conventional to index the lattice struc- 
ture by two integers m and n which define a superlattice 
translation vector Cmn = mTi + 71X2- The transverse 
momenta on an (m, n) tube are quantized to the val- 
ues 6n = So + 27rN/\Cmn\ where Sq = 27rsgn(i^)/3|Cm„|, 
depending on the sign of the chiral index of the tube, 
ly = mod(rt — to, 3). In Eq. Sx — X6n. 




FIG. 1. Planar BN forms a honeycomb lattice with B 
and N occupying alternate sites (upper left). A BN tube is 
formed by wrapping the sheet through the translation vector 
Cmn (lower left) quantizing the transverse crystal momenta. 
The upper right hand panel shows two representative lines 
of allowed momenta for a tube with a nonzero chiral index 
V = mod(m — n, 3) and chiral angle 6. The lower right panel 
shows the structure for the nonchiral armchair (10,10) wrap- 
ping of the BN sheet. 



The eigenvectors of H in Eq. (^ are 



~i4>/2 



V2 



V2 



(2) 



where E = y'v^k^ + P^) + , cj) = tun-^Sx/ \k), 
u = ^y{E - A)/E, and v ^ y/{E + A)/E. In the rest 
of this paper we will set the Fermi velocity tii;- = 1 so 
that energies and momenta are measured in the same 
units. Eq. (||) explicitly includes a gauge function 
6mk since the overall phase of the Bloch function is 
not fixed, and our calculation will require differentia- 
tion with respect to k. The Bloch eigenf unctions of our 
problem are expressed as the product of three factors 

We extend the model Hamiltonian in Eq.(|^) to in- 
clude coupling of electrons to the oscillating optical fields 
E(t) = E e"*'^* + c.c. through the dipole operator [|J|] 



Hint = — e E(<) • r. For a field polarized along the tube 
direction the matrix elements of the dipole operator be- 
tween band eigenstates m and m' at crystal momenta k 
and k' are expressed using a formulation due to Blount 

rm'm{k\k) = (?/'fe'm'k|'0fem) = {'4'k' m'\dk\i'k,n) 
+ ^S{k-k')Y^ Oafs Qm'P 

{k')d, (k) . (3) 

a/3 

The second term on the r.h.s of (||) forms the connection 
ig*n'{k) ■ dkOrnik) duc to the k dependence of the 
eigenstates in (||) |]l7| . 

In general the shift current J can be expressed in terms 
of the band off-diagonal matrix elements of the velocity 
operator v^^ (v = r]), and the band off-diagonal 
term in the density matrix Pap^ calculated to second or- 
der in the exciting fields E, as follows 

/dk 
^ V„„(fc) p„„(fc) = Je + + Jr ■ (4) 



We focus on the component of the excitation current Jo 
along the tube direction. After evaluating the sum in (Q), 
we arrive at the transparent result 

/dk ■ 
— fcv{k) Rev (k) ^ e n TZ . (5) 

Here fcv{k) is the transition rate at wave vector fc, and 
the shift vector Rcv{k) is given by 



Rcv{k) = dk9yc{k) + £,cc{k) - Cvvik) 



(6) 



with ^„c(fc) = \^vc{k)\ e'^-C'). The shift vector Rcv{k) 
is invariant under the gauge transformations exp(i0mfc) 
and we will evaluate it in a gauge with dkOmk = 0. 

For an incident electric field polarized along the tube 
direction, interband excitations are allowed only between 
band states with the same transverse momentum S. We 
consider transitions between the lowest two azimuthal 
subbands {N — 0). Then using the eigenfunctions in 
Eq.(|) we find 



£,vc{k) = 



A 



rdkE + idk 



^£:2 - A2 

which gives the off-diagonal contribution in Eq. (^ 



dk&vc = 9/0 tan" 



kA 



The diagonal elements in Eq. (J^ are 



^ A 
Okgm —^2E 



S'^ + k^ 



(7) 
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(9) 
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where the sign is negative for valence states and positive 
for conduction states. Combining Eq. (^) and Eq. ^ 
we obtain the shift vector 



Rcvik) 



25A 



(,52 + fc2)VFTA2~rP 



(10) 



Comparison with Eq. (|^) shows that the shift current 
(Eq.(^) opposes the direction of the ground state polar- 
ization of the tube. 

This shift vector is odd in the symmetry breaking po- 
tential A and odd in the transverse momentum S. There- 
fore, two BN tubes with nearly the same wrapped lattice 
structure but opposite chiral indices v exhibit reversed 
ground state polarizations and opposite shift currents. 
Armchair BN nanotubes with wrapping indices (m, m) 
have chiral index 1^ = and do not exhibit a longitudi- 
nal shift current (as shown in the middle panel in Figure 
^ for the (5,5) tube). Zigzag BN tubes with wrapping 
indices (m, 0) can be grouped into three families with 
positive, negative or zero shift currents distinguished by 
the sign of the chiral index i/ (an example for the (8,0) 
tube with = — 1 is in the left panel of Figure |^) . It is 
important to notice that an isolated flat BN sheet has a 
threefold symmetry axis perpendicular to the BN plane 
and therefore it has zero electric polarization in its tan- 
gent plane, by symmetry. Thus it has no shift current for 
excitations at the band edge. The nonzero shift vector 
in Eq. (10) when A: ^ is a remarkable long wavelength 



quantum mechanical effect that reveals the quantization 
of the transverse momentum on the wrapped tube. 




(8,0) 



(5,5) 



(7,3) 



FIG. 2. Electronic shift currents produced by excitation 
of BN nanotubes with light polarized along their axis. For 
the (8,0) wrapping the chiral index is — —1, and the shift 
current is purely longitudinal. The (5,5) armchair wrapping 
has chiral index v — and has zero longitudinal current and a 
nonzero azimuthal current which runs "counterclockwise" on 
the tube surface. The (7,3) wrapping has a chiral shift current 
with nonzero longitudinal and azimuthal components. 

The shift current Jo in Eq. (||) is determined by the 
average shift vector TZ, describing the change of polarity 
in the excitation, and by the carrier injection rate h. The 
BN bond has significant ionic character with |]l8| 2A « 4 
eV. This gives for excitations near the band edge, fc ^ 0, 
a shift vector with a magnitude TZ w 0.06 nm for a (17,0) 
tube. For an incident intensity S = 100 kW/cm^ we can 



also obtain the interband excitation rate; we estimate 
that n « 70 nm~^fis~^ produces a net electronic shift 
current of Jc ~ 0.67 pA. This value is slighly larger than 
one obtains for excited bulk polar semiconductors. In 
a coaxial multiwall tube Rev can change magnitude and 
sign on each wall. Although opposing currents on differ- 
ent tube walls reduce the macroscopic current a complete 
cancellation is very unlikely since the magnitude of Rev 
depends on the tube radius. We also expect that the re- 
combination current Jr will oppose Jo since both currents 
rely on optical transitions with longitudinal polarization. 
However, this should not prevent the observation of Je 
using pulsed laser excitation since the timescale for Jj. is 
of order nanoseconds. 

This formalism can be extended to study chiral shift 
currents with nonvanishing longitudinal and azimuthal 
components. Interestingly, the long wavelength Hamilto- 
nian in Eq. (|l]) gives zero azimuthal current for any chiral 
index 1/. As Q± = k zL iS ^ this occurs because of a 
cancellation of contributions from the K and K' points 
whose long wavelength Hamiltonians have opposite hand- 
edness. However, azimuthal currents arise due to discrete 
lattice effects occurring at higher orders in an expansion 
in Q±. To study this, we define two unit vectors at each 
point in the tube surface, perpendicular to the local out- 
ward unit normal n: e^ which points along the tube axis, 
and Cy = n X Cx which points "counterclockwise" on the 
tube surface. A chiral angle 9 is defined as the angle 
between and the Ti lattice direction as shown in Fig- 
ure |l| (with this convention armchair nanotubes have 9 
— 0). Then the lowest nonvanishing contribution to an 
azimuthal (y-polarized) shift current must have the form 
cos(6</) — 36') where = tan~^((5/fc). The dependence 
on 60 means that azimuthal currents first occur in the 
continuum theory at sixth order in \Q\. 

We can study these terms by replacing the kinetic 
energy terms in the Hamiltonian of Eq. (^ by the 
discrete lattice counterpart ivpik + iS) — ivFQ+ — *■ 
t^^e^'^-^^'^''"^" where t is the hopping amplitude be- 
tween sites connected by the bond vectors r„ . A system- 
atic expansion of the azimuthal shift current in powers 
of Q± then gives contributions with the correct lattice 
symmetry and shows that the ratio of the azimuthal and 
longitudinal in-plane current densities scales as 



Jy 



(fc2 



^ Jx <5(A2 + fc2+^2) 



COS 



•-39). 



(11) 



The azimuthal current for any zigzag wrapping is zero 
because of the angular factor, while its longitudinal cur- 
rent, shown in Figure^, is determined by the chiral index 
v. For armchair tubes and all non-zigzag tubes with chi- 
ral index v = 0, we have S — 0, so the shift current 
is purely azimuthal with a direction determined by the 
sign of the potential A. These energy dependent currents 
vanish at the band edge as oc (i?2 — A2)'^. Non-zigzag 
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tubes with chiral index v ^0 have a nonzero 5 and thus 
the azimuthal current is nonzero with p cx (5/A)^((5t)^. 
There the shift current circulates in a helical pattern on 
the tube surface as shown in the right panel of Figure 
^, with a chirality determined by the index v. The 
pitch is also energy dependent and increases according 
to p{E) = p{0) + a{E'^ - A^). Thus, as shown in Fig- 
ure ^ the zigzag nanotube generate a photocurrent like a 
wire, the armchair tubes like a coil, and the chiral tubes 
can exhibit a helical current on the tube surface. 

Alternatively, the shift current can be studied in a 
semiclassical model, by considering an expansion in t/ A, 
where t is the nearest neighbor hopping amplitude, and A 
is the site diagonal potential. In the limit t/ A << 1 the 
tube can be regarded as a network of independent bonds 
whose excitations are superimposed. Here the local sym- 
metry of the zigzag tube requires zero net azimuthal 
current while the armchair wrapping has a nonzero az- 
imuthal current with its direction determined by the the 
sign of A. However, it is difficult to correctly describe 
the physics of the longitudinal current using this model. 
This is because the longitudinal current is controlled by 
the quantization of the transverse momentum S. The first 
terms in the semiclassical expansion that are sensitive to 
S occur at order (t/A)^ where N is the number of bonds 
(« 40) around the tube circumference. This reflects the 
fact that the longitudinal current is fundamentally a non- 
classical quantity for this system. 

The longitudinal shift current reduces the ground state 
electric polarization of the nanotube. This can lead to a 
fast readjustment of atomic positions on the tube walls. 
Generation of voltages by a mechanical elongation 
of NCS nanotubes is a related effect. In both situations 
the response can be stronger in buckled nanotubes pC| ] 
or if different ty pes of atoms occur in different layers, as 
in M0S2, WS2 in or GaSe nanotubes H. Recently, 
mechanical motions, strikingly similar to those found in 
some biomolecular systems |23| ], were observed in irradi- 
ated nanotube bundles Interesting mechanical re- 
sponse to an applied dc bias, is observed in carbon nan- 
otubes immersed in a solution of NaCl |^^. Here, the 
motion is caused by attraction of ions of a given polarity 
to the nanotube. The speed but short duration of the 
shift current-induced mechanical response complements 
the slow but steady state effect of a dc bias. 

Experimental observation of the shift current in BN 
nanotubes would provide a striking illustration of funda- 
mental concepts in the modern quantum theory of polar- 
ization. The effect is also promising for new applications 
since it couples the microscopic physics on short length 
scales, which can be tuned by local fields and mechanical 
loads, to the long distance properties of the system (po- 
larization, photovoltage, etc.) Finally, there is a rapidly 
growing family of related submicron one dimensional ma- 
terials to which these ideas can be applied. 
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